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Based on the two Higgs doublet model, we study the effect of Higgs- 
boson exchange on the (super)heavy quarkonium QQ, which induces a strong 



£SJ ■ attractive force between a (super)heavy quark Q and an antiquark Q. An 

interesting application is the decay of (super)heavy quarkonia QQ into a Higgs 
boson associated with gauge bosons. The criterion for making the QQ bound 
state is studied. We also show that non-perturbative effects due to gluonic 



field fluctuations are rather small in such a heavy quark sector. Possible 



enhancement for productions and decays of QQ bound states made from the 



fourth generation quark Q is discussed for pp (at the Tevatron) and pp (at 
the LHC) collisions. 
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I. INTRODUCTION 



A study on quarkonia T(QQ) composed of a (super)heavy quark Q and an antiquark 
Q (a possible and interesting candidate of Q is the up(£/)- and/or down(D)-quarks in the 
fourth generation family) is required in current particle physics for testing the standard 
model (SM) and/or searching for signals for physics beyond the SM. It is nowadays one of 
the most interesting subjects since the subject can be studied, with high priority, in the 
forthcoming experiments at high energy hadron colliders, i.e. the Tevatron and LHC. In 
particular, no theoretical arguments are seen to rule out the (super)heavy quarks and the 
(super)heavy quarkonium states with the masses around hundred GeV or even a few TeV 
(see, e.g. [1]). As the recent example, in one of the extended models, the little Higgs model 
[2], there could be an additional heavy quark q with a mass of order 0(1 TeV) to promote 
the quark doublet q to quark triplet under the global SU(3), Q = (q,q). By preserving 
the global symmetry of the coupling, the one-loop quadratic divergence to the top quark is 
removed. In fact, the SM and its extensions, e.g., the Minimal Supersymmetric Standard 
Model (MSSM), do not explain the family (generation) structure of the quark masses. Each 
quark has an arbitrary Yukawa coupling and hence is independent of the family to which it 
belongs. It is required to explain the family structure and the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix for the quark sector in any extension of the SM or even in the SM. 
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It is known from the history of particle physics that the first signals for c- and 6-quarks 
in hadronic collisions were leptonic decays of their J/ip{cc) or T(bb) bound states. Can (su- 
perheavy quarks also be first discovered through the decay of their quark-antiquark bound 
states into lepton pairs? It seems the answer is apparently not transparent, because one 
of the main properties of (super)heavy quarkonia is concerned with the appearance of new 
decay modes into weak bosons and even Higgs bosons in the final states. Once the Higgs 
boson H is discovered, one needs to measure its couplings to other particles. The value 
of i7-boson couplings can be extracted by measuring a variety of Higgs boson productions 
and their decay modes. Thus, it is important to find the if-boson in as many channels 
as possible, including its production coming from the decay of (super) heavy quarkonia. 
Such a Higgs boson can be looked for in the following new decay modes; T{QQ) — > HZ, 



T(QiQi) -> HT{Q 2 Q 2 ), T{QQ) -> 7 i7 and T{QQ) -> ggH^H, where T{QQ) carries 



quantum numbers of J PC = 1 , 0~ + . Production of heavy quarkonia like bb and it associ- 
ated with a Higgs boson emission in the decay of extra gauge bosons Z' has been studied 
in [3,4]. Here we are interested in the case in which the new decay modes mentioned above 
become dominant and hence the branching ratio for a single heavy quark decay accompa- 
nying the real weak boson emission Q — > qW (q is a lighter quark), leaving q as a spectator, 
is small. Of course, it is necessary to examine whether the spectator mode can be dominant 



The direct decay of the top quark, t — > Wb [5,6], and the flavor changing top quark de- 
cays, t — > erf {eg, cZ) and t — * cH [7], in the both SM and two- Higgs doublet model (2HDM), 
have been studied intensively for the last decade. The SM predictions of the branching ratios 
{BR) for those decays t — > c~y, t — > eg, t — > cZ and t — > cH are significantly small, being 
BR ~ 5 x KT 13 , 4 x KT 11 , 1.3 x KT 13 and 10 -14 - KT 13 , respectively. Obviously, the 
rare decays are out of interest here since they are very difficult to be observed in hadron 
colliders even at the highest luminosity and thus we neglect them in this work. The sensi- 
tivity of current experiments at Tevatron Run II or future experiments at forthcoming LHC 
closely approaches the rate required for ruling out the Higgs boson production or discover- 
ing it through the decay of (super) heavy quarkonia mentioned above. In this connection, 
precise theoretical estimates of the rates are required for an unambiguous interpretation of 
experimental upper limits. 

One cannot exclude the possibility of the new strong interactions which primarily control 
the dynamics of (super)heavy quarks such as the 4th generation up(£7) and/or down(TJ) 
quarks. In one of the "top-color" models [8] with 1 TeV scale, there is the following "top- 
color" gauge structure 



where SU{3) A x U{l) n , SU{3) h x U{l) Yh and SU{3\ x U(1) Y , generally coupled to the 
4th, 3rd and first two generations, respectively. The U{l) Yi are just rescaled versions of 
electroweak U{1) Y into the strongly interacting world. In this model, below the symmetry- 
breaking scale fJ>sB, the spectrum includes massive "top-gluons" , which mediate vectorial 
color-octet interactions among (super)heavy quarks Q{= t, U, D) 



or not. 



SU{3) 4 x SU{3) h x SU{3) t x U{l) n x U{l) Yh x U{l) Yl x SU{2) L 

-> SU{3) QCD x 1/(1) 

EM j 



(1) 




(2) 



2 



If the coupling k lies above some critical value n cr it, the heavy quark condensate (QQ) can 
be formed. The strong "top-color" dynamics can bind Q and Q into a set of "heavy-pions" 
(QQ). The criterion for existence of a (super)heavy quarkonium is that the binding en- 
ergy tB should be larger at least than the total decay width T tot of its quarkonium, namely 
c = (r tot /es) < 1. Since such a (super)heavy quark-antiquark bound state is considered to 
be a non-relativistic system, the quark potential model should be applicable to the analysis. 
Then as pointed out in [9-11], one cannot exclude the possibility of the Higgs-boson interac- 
tion which dominates significantly over the one-gluon exchange ~ -(4/3)a s (mg)/r for the 
(super)heavy quarkonium, where a s is the strong coupling constant depending on the quark 
mass ttiq and r is a distance between a quark and an antiquark. We show in Sec. 2 that 
the strong binding force due to the Higgs-boson exchange gives rise to a necessary condition 
c < 1 for enabling the heavy quarkonium to exist and leading to observation of its resonance. 

On the other hand, in the physics of interplay among quarks, it is well known that the 
exact QCD vacuum should contain the fluctuations of gluonic fields at large scales [12]. 
Those non-perturbative fluctuations cause the distortion of interactions between quarks and 
antiquarks. We consider (super)heavy quarks as external objects allocated in the gluonic 
vacuum. Here, we study those non-perturbative fluctuation effects of the gluon field on the 
decay of a quarkonium T(UU) into the Higgs- and Z-bosons. Our result is based on the well- 
known statement (see, e.g., [12]) that the non-perturbative effect on dynamics of heavy quark 
systems is expressed in terms of vacuum expectation values of the local operators constructed 
from gluonic field operators. The leading effect is proportional to a matrix element of the 
form (GV 2 (0)) , where G°„ (x) is the standard gluonic field strength tensor with color 
indices a = 1, 2, • • • , 8. The lowest level of the (super)heavy quarkonium is determined by 
the color-singlet Yukawa-type attractive force mediated by the "light" scalar x-boson. Here, 
we are interested in the corrections due to non-perturbative gluonic fluctuations in the exact 
QCD vacuum. In this paper, we show that in a (super)heavy quarkonium decay such as 
T(UU) — > hZ (h means a lightest CP-even Higgs-boson in 2HDM), the non-perturbative 
fluctuation effect of the gluonic field can be calculated, to some extent, without detailed 
knowledge of the vacuum structure and furthermore, it gives a negligible result. 

The outline of this work is organized as follows. In Sec. 2, we discuss an effective 
potential mediated by a Higgs boson. An effective model for the lower-energy theorem will 
be discussed in Sec. 3. Finally, in Sec. 4, we give our conclusion and discussion. 

II. EFFECTIVE POTENTIAL VIA HIGGS BOSON EXCHANGE 

First of all, one should notice that the total cross-section for the process pp(pp) — > QQ 
is strongly related to the underlying subprocesses. For example, at =14 TeV, the cross 
section for various processes of Q and Q pair productions with Q = U and/or D in the mass 
region my ~ mo — n^4 = 0.2 — 0.6 TeV (m^ is the mass of the fourth generation quark), 
becomes a ~ 1CT 1 — 1CT 2 (pb) and a ~ 1 — 10~ 2 (pb) for 7/Z production and W production, 
respectively, in the q and q annihilation channels, and a ~ 10 — 10 _1 (pb) for the two-gluon 
fusion channel. 

Let us consider a production (mpp orpp) of (super)heavy quarkonium T(QQ) followed by 
the decay process T(QQ) — > hZ with Q = t,U, D, being assumed to be the dominant decay 
process of T(QQ). Then, it is supposed that this dominant mechanism at high transverse 
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momentum involves, e.g. the production of a gluon that produces a color-octet Q and Q pair 
which then fragments into a color-singlet bound state T(QQ) by emitting two or more soft 
gluons. This T(QQ)-state should be transversely polarized at high transverse momentum, 
since it is emanated from a gluon which has only transverse polarization states. 

In the lowest bound state QQ, the quark Q and the antiquark Q are assumed to be 
located at a distance 

r ~ [m Q A(m Q )] _1 (3) 

which is small compared to the scale of strong interactions. A(mg) is a strength of the 
interaction between a quark and an antiquark. The wavefunction of the lowest bound state 
is proportional to exp(— /ir) with /i ~ mg A(m(j). Note that for r = fi' 1 being smaller 
than in the typical size of fluctuations, our approach becomes applicable since the potential 
dumps exponentially for a distances r ^> 

Let us consider a simple model (Model I) where the dominant effective potential for a Q 
and Q system at a small distance looks like [10,11] 

V eff (r) ~ a s {m Q ) - ^^o) exp(-m x r) , (4) 

with 

2 

X ( m Q^xQ) = {^e xQ , (5) 

and £ x q reflects the model "flavor" in the strength of the interaction between the scalar 
X-boson and a heavy quark Q (£ x q = 1 in the minimal SM, otherwise £ x q > 1). v is 

the vacuum expectation value of Higgs boson, v— 246 GeV (in 2HDM, v = \Jv\ + t|, v i 
and t>2 are two neutral Higgs field vacuum expectation values) and Cf is the color factor, 
Cf = 4/3 for color SU(3) group. In Fig.l, we show the ratio of the combined coupling 
otcomb = (4/3)a s (m Q ) + \(m Q ,£ xQ ) ■ exp(-m x r) (see (4)) to the pure QCD coupling a s = 
(4/3)o; s , as a function of a heavy quark mass tjiq for different values of £ x q and m x . The 
second term in a com b is appropriate for r ~ [mgA^g)]" 1 . The ratio becomes somewhat 
bigger for smaller Higgs-boson masses. 



Because of our demand, A(mg,^Q X ) > CVa: s (mQ), for a relevance of the x-boson inter- 
action, the lower bound on uiq is given as 



m Q > -^(AnC F a s yl\ (6) 



which leads to tuq > m t even if £ x q= 2 (see Fig. 2). 
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a) m Q [Tev] m Q [Tev] 

Fig. 1 Ratio of the combined coupling a com b to the pure QCD coupling a s as a function of a 
heavy quark mass for (a) £ x q = 1 and (b) £ x q = 2, 3. The curves are presented for x-boson 
masses m x = 90, 100 and 120 GeV. 



The requirement of the positivity of the variational parameter 



^^l AW f?l (7) 



2 (Am Q ) 2 + 2 



x 



entered in both the bound state wavefunction \l/(r) = 2 // 3 / 2 exp(— /xr) and the binding 
energy (see, for details, [11]), 

eg = 2A /i ; (2/i -^ ) (8) 
leads to an upper limit on m x (see Fig. 3). 



In the Model I, the ratio c = (T tot /e B ) < 1 could be guaranteed for T(UU) quarkonia to 
be formed by a strong attractive force via scalar Higgs-boson exchange with a sufficiently 
"hard" Yukawa coupling A(mQ,£ x g). The total decay width T tot is given by a sum of two 
terms 

T tot = T T + T U , (9) 
where the width Ft is defined by the following decay channels 

T(UU) -> hZ, -yZ, jh, W + W~,bb, it, T + T-,n + fT, ggg , 
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Fig. 2 Lower bound on (super)heavy quark masses as a function of £ x q. 

and the single quark decay width IV is a sum of the following contributions: U — > 
DW + ,bW + ,bH + {H+ is the charged Higgs-boson in 2HDM). In the decays of T ([/[/) pre- 
sented above, the radiative channel is suppressed by the coupling constant a, the ggg channel 
gives the small contribution due to factor, and furthermore, the productions of the pairs of 
quarks and antiquarks or leptons and antileptons are also small because they follow via two- 
loop diagrams, where the amplitude is very small because of the presence of a (intermediate 
photons) or the Fermi constant G F (virtual VF ± -bosons). 

We do not consider the contributions from the decay U — > bH + because it is expected to 
have a rather small probability to be observed as expected from the following consideration: 
the Higgs-boson mass sum rule 



with one-loop correction 5 < 10% [13] does not allow the production of charged Higgs- 
bosons in the decay of top quark and, perhaps, also of U quark due to the kinematical 
reason in the decoupling limit, (m^/m^) <C 1 (m A is the CP-odd Higgs-boson in 2HDM). 
In addition, the experimental data at the Tevatron do not yet clarify the status of t — > H + b 
decay. The CDF results in direct search of a r-lepton emission from top quark decays give 
an upper limit on the branching ratio BR(t — > H + b) ~ 0.5 — 0.6 at 95% C.L. in the range 
60 GeV < m H + < 160 GeV, assuming BR(H + -> tv t ) = 1 [14]. Furthermore, the DO 
Collaboration excludes BR(t -> H + b) > 0.36 at 95% C.L. in the region 0.3 < tan/5 < 150 
and m H + < 160 GeV [15]. Assuming J2x=h+,w+ BR(t — > bX) = 1, the decay width of 
the if + -boson channel in U- and t-quark decays is small and hence this channel is out of 
interest. 

Therefore, for the case of U quarks, one can expect that the T tot is given by 



m 2 H+ ~ {m\ + m 2 w ) (1 + 8) 



(10) 



T T (T(UU) -> hZ, W + W~) + YuiU -> DW + , bW + ) . 



(11) 
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Fig. 3 Upper limit on the scalar x-boson mass as a function of tuq for different values of £ x q. 
The regions above the corresponding curves are excluded. 



The main contribution to IV arises from the channel T(UU) — > hZ, whose decay width for 
T(UU)(1 ) is given by the following expression (see, e.g., [11]) 



T{T{UU)^hZ) = ^r 1 l u 



x 2 , 1 2 (™U 

2 \m w 



with 



aw = en/sly ,a z = aw/c w (sw = sm9 w , cw — cos 9w), 



r = (l-g^r) • <I» = ( I 



a 9 \ 1/2 

4m|\ 
M| J 



f = /(A,m [/ ,m x ) 



(Amu) 2 - m 2 x 
(Xmu) 2 + 2ml 



, Mt — 2 . 



(12) 



mjy is the mass of an up(£/)-quark of the fourth generation. We suppose that the couplings 
of the Higgs-boson h and the £/-quark have the same form as those for the couplings of the 
Higgs-boson h and the top-quark (see the review [16]) 



m 



r] hU ~ H J" s 2/3C 2/3 tan /3 , 



(13) 



where 82/3(02/3) = sin2/9(cos2/3) and tan/5 is the standard ratio between two vacuum ex- 
pectation values for two Higgs doublets in 2HDM. 

Using the standard formulae [5,6], the single £/-quark decays are given for U — > DW + 
channel by 



r u (u^DW + ) = ^^-\v Ii 12 
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and for U bW + decay by 

Tu(U -> &W + ) = r (c/ 

where 



bW + ) (1 - A) , 



r (£/" -> w + ) 



G F mfj 



V Ub \ 2 fy(3-2pt v ) 



(14) 



(15) 



(16) 



with p w = \Jl — m^/mfj. Vud and Vub are generalized CKM matrix elements and the 
0(a s ) QCD correction [6] in (15) is 



A = 



C F a s 
2tt 



2tt 2 



(17) 



As an example, for a typical set of parameters: m\j= 400 GeV, mu= 300 GeV, m x ~ m/j= 
100 GeV, |V{/£)| ~ 1, ^ x [/ = 2 (?7^[/ — 1) we give the numerical results for those decay widths 

I2rj 2 hu GeV, 

(18) 



r(T(UU)^hZ)~2.2 V 2 hU GeV, 

T(U -> ,D^ + ) ~ 1.21 GeV, r(C/ -> W + ) ~ 0.0021 GeV, 



where c = V tot je B = 0.27 at the calculated value e# = 12.65 GeV. 

The contribution of the decay U — > is negligible due to rather small value of 

\V Ub \ ~ 10~ 2 . We have checked the possibility of the existence of bound states T(QQ) 
composed of U (and -D)-quarks, starting at the lowest value of £ x q > 1, where the increasing 
£ x q gives rise to an effect on c" 1 > 1 (see Fig.4). 



The expected event topology of the decays U — > bW + and U — > DW + are similar to that of 
t — > bW + . However, for the down-type heavy quark D, the dominant decay mode could be 
D — > tW~ and thus, in the case of Z)D-pair production, the final state can consist of two 
pairs of leptons and neutrinos (originated from decays of the t quark and W boson) with 
different flavors, in general. Now, by taking, as an example, the following parameters: m£>= 
400 GeV, m x ~ m h = 100 GeV, £ x d = 2, \Vot\ — 0.012, we can obtain the following decay 
widths 

T(T(DD) -> hZ) ~ 2.25 rg D GeV , T(D -> tW) ~ 1.47 MeV . (19) 

One can see that the decay width T(T(DD) — > /iZ) can be enhanced by the Yukawa-flavor 
factor (within 2HDM) rjh d = — sin a/ cos/? at large values of tan/5. 

For comparison, let us consider an instructive example, i.e. the decay of T(UU)(0~ + ) 
states to h- and Z-boson, where the decay width is given by the following expression (see, 
e.g., [11]) 
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Fig. 4 Ratio c = T tot /e B being responsible for occurrence of the bound state T{UU) as a 
function of the quark mass my at different values of (, x u for m x = 100 GeV. The region 
below the value c = 1 is allowed. 



V(T(UU) _ hZ) = ^^„<4 g^f (20) 

The numerical estimation gives T(T(UU) — > /iZ) ~ 6.80 t^c/ GeV, which is roughly 3 times 
larger than that in the case of T(UU)(1 ) — > hZ decay mode given in (18) and hence leads 
to larger production rate of T(UU)(0~ + ) than T(UU)(\ ). This simple example confirms 
our belief that the most promising candidate of the (super)heavy quarkonium which could 
be searched at the LHC should be the pseudoscalar state T(UU)(0~ + ). 
It is interesting to estimate the effect of the scalar x-boson exchange on the (QQ) production 
cross-section at different m x as a function of itlq. 

If the Q-quark is relatively long-lived, the peak of the cross-section at a QQ resonance 
due to the one-gluon exchange is given by 

'-*(%)-. (21 » 

for a given center-of-mass energy s. The cross-section (21) has a strong sensitivity to a s and 
decreases sharply with increasing itlq because of the rapid grougth of Tq (see Eqs. (14) and 
(15)). To leading order, the scalar x-boson exchange effect is taken into account simply by 
making the replacement 

a s -> a s + a(m Q , m x , £ xQ ) , (22) 
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where 



3 mj 



16 71 V 



£v Q exp(-e) 



2 >X 



(23) 



with e = m x /{rriQ A) < 1. Here, we see a simple increase of the coupling strength between 
Q and Q. We find a relative enhancement of the effective cross-section cr e // 



-i ~Q O ~ o o ^ 

1 + 3 a « s + 3a a s + a a s 



(24) 



due to the x-boson exchange effect. In Fig. 5 we show the ratio of the effective cross-section 
(24) to the cross-section (21) as a function of itiq. The curve in Fig. 6 is appropriate for a 
massless x~ boson. 
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Fig. 5 Ratio of the effective cross-section a e ff (24) and the cross-section a c at the IS resonance 
QQ as a function of toq. The curves are for £ x q= 1.05, 1.10, 1.15 and m x = 100 GeV at 
a s = 0.118. 



In the end of this Section we give briefly the results of the decays $ — * T(QQ) + 7 
where a Higgs-boson $ ($ = h or $ = H) with 4-momentum and the mass decays 
in the heavy vector (1 ) quarkonium T(QQ) (Q = 6-quarks, t— quarks, ...), carrying the 
momentum = 2p^ (P 2 = m|), and a photon with the momentum squared k 2 = 0. 
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0.1 0.2 0.3 0.4 0.5 



m Q [Tev] 

Fig. 6 Ratio of the effective cross-section a e ff (24) and the cross-section a c at the IS resonance 
QQ as a function of tuq. The curve is for the massless x-boson. 



The amplitude of the transition $ — > T7 is [17] 

A($^T 1 ) = -2V^g v e Q ^v*Q 1 _^ /mQ y £ a k (3 (P*fa-P^ a ), (25) 

where cq is the charge of the quark Q, being the polarization vector of T(QQ) and gy is 
defined in the standard manner 

(T(QQ)\Q^Q\0) = m 2 T g v( f )f , 

and can be estimated from the leptonic decay width T(QQ) — > II (I — e , /i , r): 

4 

T(T(QQ) -> ZZ) = - 7r (a e Q 5fy) 2 m r . 

At an arbitrary large values of y = {m<$,/2rriQ) 2 > 1 there are the corrections to the amplitude 
(25) due to the one- loop 0(a s ) gluon contributions 

l-F(y)C F a a (m%-4m 2 Q ), 

where the function F > was calculated in Ref. [17]. For the decay processes h — > T(bb) 7, 
H — > T(bb) 7 and — > T(tt) 7 we are considered here as the promising channels in a wide 
range of the Higgs-boson mass we expect the relative decay width compared to final state 
quarks QQ, which is given by [17] 
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f = BR(<$> 



t(qq) 7 ) r($^T(gg) 7 ) 



BR($ -> QQ) 

647rQ 4 2 ( m Q\ 2 tSL, n 

^^£y K {y,a s ) 

3 \m$/ 



r($ - qq) 



1-1/2 



(26) 



where 



K(y,a s ) ~ 1 



Q! s (m 



7T 



C F In 2 ln(4y) 



will be important for the processes h — > T(66) 7, iJ — > T(fefe) 7, where (m(,/m$) 2 — > 0. 
In Figs. 7 and 8 we plot the relative widths of /i — > T(66)7 and — > Y(66)7 decays, 
respectively, compared to 66 final state, vs. and mjj, respectively, and if — > T(tt) 7 
decay (see Fig. 9) compared to tt-state as a function of win- 

The recognizing of the quark- ant iquark bound state can be shown through the resonance 
structure having a specific signal, e.g., the final leptonic pairs e + e~, and t + t~. This 

resonance is expected can give one clear evidence of a bound state production over the QCD 
background if the production would be significantly large. 
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Fig. 7 The relative decay width T(h 
to h — > 66 channel. 
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Fig. 8 The relative decay width T(H — > Tibb)^) as a function of iJ-boson mass m H compared 
to H — > bb channel. 

III. EFFECTIVE MODEL FOR LOWER-ENERGY THEOREM 

Before proceeding to the lower-energy theorem applied to the transition T — > hZ (Model 
II), we give the general expression for the amplitude of the decay mentioned above [18] 



A(T -> hZ) 
1 

v 



<z\£int\<t>T) 
Q=t,U,D 



12tt 



Vh Q G%G a ^- J2 Vh q m q qq 

q=u,d,s,c,b 



(27) 



Here, \4>z) and \4>t) are eigenstates responsible for Z-boson and (super)heavy quarkonium 
T = T(QQ), respectively. The second term in (27) corresponding to the contribution of 
light quarks comes directly from the interaction Lagrangian 



J2VhimiU+ J2 Vh q m q qq+ ^ Vhqm Q QQ- 

q=u,d,s,c Q=b,t,U,D 



-2 VhW m 2 w W+W>*- 



2 ry2 

r) hZ m z Z 



hi > 



(28) 



while the gluonic contribution given by the operator in the first term of (27) arises from the 
coupling of the Higgs boson h to t-, U- and D-quarks through the standard loop mechanism 
(see, e.g., references in [18]) with Nh number of heavy quarks in the loop. The functions 
rj hi in (27) and (28) (i=leptons (/), light quarks (q), Q, W ± , Z) are model-dependent ones 
describing a deviation from the SM picture where all of r\ hi are equal to unity. Obviously, 
considering the fermionic sector, only the terms containing quarks with masses mq > m h 
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f(H ^T(tt)y)x10 3 




m H [Tev] 

Fig. 9 The relative decay width T(H — > T(it)~f) as a function of f/-boson mass m H compared 
to H — > it channel for different values of g v . 



are relevant for giving rise to the amplitude provided by the matrix element of (27) where 
the only gluonic part survives. Our aim is to calculate the amplitude (27). Fortunately, 
QCD gives the trace of the energy-momentum tensor 0^ in the following form [18]: 



G%CT* + '£ mq qq + 'E mQ QQ, 
Q Q 



(29) 



where b$ is the first coefficient of the /3(a s ) function in QCD with 6 flavors of quarks. At 
the zero momentum transfer, the matrix element of (29) between any different eigenstates 
X and Y of the Hamiltonian is vanishing [19,18] 



(x\e^)\Y) = o 



(30) 



with (X\Q 00 (qQ,0)\Y) = 0. Using the condition (30), one can obtain the following relation 
between the matrix elements containing the gluonic part and the heavy quark terms: 



"6 a s na n n V a 



t) — \4>z mb bb + m t tt <fi T 



(31) 



where the light quark contributions are neglected. Taking into account the gluonic anomaly 
effect through the replacement of the top-quark mass term [19] 



2 a 

m t it^---^G^G^ a , 
6 on 



(32) 
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Eq.(31) transforms into 



^5 a s na siiiva 



>T = \(PZ 



nib bb 



(33) 



where 65 = b% + 2/3 . Let us consider the amplitude (27) in the following form 



A(T -> hZ) 



1 

v 



E 

Q=t,U,D 



lYvr 



(34) 



where the terms containing u, d, s and c quarks were omitted because of its small mass effect 
on the intermediate quark loop. Comparing Eqs. (34) and (33), we find the amplitude for 
the decay T — > hZ in the limit of vanishing of the four-momentum of the Higgs-boson h: 



K 


4>z 


Us 1 






8tt \ 



A(T -> hZ) = — 



where the couplings r] h q are [16]: 

cos a 

Vht(U) 



sin/5 



° Q=t,U,D 



sin(/3 — a) + cot /3 cos(/3 — a) , 



(35) 



(36) 



VhD 



sin ct 

cos /3 



sin(/3 — a) — tan /3 cos(/3 — «) 



(37) 



For numerical estimation, we use the decoupling limit where Eqs. (36) and (37) are trans- 
formed into the following ones: 



Vht(u) ^ 1 + z s 2/ 3 c 2/3 tan 1 f3 , 
VhD ^ 1 - z s 2/ 3 c 2/3 tan j3 , 



(38) 
(39) 



where z = {m z /titia) 2 ■ To calculate the amplitude (35), one has to estimate its matrix 
element as those given by the soft non-perturbative gluonic field. This contribution occurs 
as the excitation of the non-perturbative gluon condensate in an environment of a pair of a 
quark and an antiquark bound at the scale 



r- 1 
bound state 



A Q = m Q a s (fi ~ m Q a s ), 



(40) 



which is larger than the scale of strong interactions A with « S <1. This non-perturbative ef- 
fect can be estimated in the transition T — > Z within a minimal point-like source F a s E 2 (x) 
[19] with E being the electric component of the gluonic field, and an arbitrary constant F 
defines the strength of this source. Hence, the only remaining work is to calculate the 
following two-point function 



W{G au ) 



dxe iqx (0 



T{Fa s E 2 (x), 



4a 







(41) 
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embedded in the non-perturbative amplitude A^p of the decay T — > hZ 

Anp(T _ hZ) = J_ {, _ | 1 ^ , kQ ) . W(G tll/ ) . (42) 

Here, (3{a s ) = {—b 5 a 2 a /2ix) + 0(a 3 s ). Using the tricks followed by the authors in [12,20], the 
calculation of the amplitude (42) can be achieved in the framework of the QCD low-energy 
theorem approach in the limit q 2 — > 

A M T -* HZ) = ^ (l - £ ^ E ^ ^) C)„ • (43) 

where {(a s /ir) G 2 IU ) Q is the standard gluonic condensate. Using these formulas, now we can 
estimate the corrections due to non-perturbative fluctuations of the gluonic field in the decay 
T(UU) -> hZ 

T(T(UU) -> /iZ) = T (T(UU) -> hZ) (1 + <W) , (44) 

where T (T(UU) — > /iZ) is the decay width (12), while the non-perturbative correction 
factor Snp is defined as 5^p = rArp/r with 



i 2 \ 1/2 

T NP (T(UU) ^ hZ) = — — \A NP \ 2 ( 1 - ^ ] • (45) 



1 14 I 2 h 4 < 

1« ^ATP 1 2 

IburriT \ raj. 



The <5tvp correction is obtained at the order of the magnitude of 10 8 for the strength 
parameter F — 1. 



IV. CONCLUSION AND DISCUSSION 



In this paper, we have concerned ourselves with the question of existence of (super)heavy 
quarkonia and their decays with production of a Higgs-boson, e.g., T(UU) — > hZ. The 
possible existence of the (super)heavy quarkonium has been studied in the framework of a 
simple Higgs-boson potential model (Model I) having a phenomenological "hard" Yukawa 
couplings \(rriQ, £xq) between the (super)heavy quark and the Higgs-boson x- Furthermore, 
to estimate the decay width of T(UU) — > hZ by taking into account the fluctuations of 
the gluonic field (Model II), we used the conformal properties of QCD. We found that 
the effect of non-perturbative fluctuations of the gluonic field in the total decay width 
Y(T(UU) — > hZ) is rather small. T^p is quite sensitive to the strength F of the gluonic 
point-like source, since it is proportional to square of this strength. The constant F can be 
estimated within the approach of the Model II through the decay, e.g., T(UU) — > II (I — fi, r) 
with A NP ~ (ll\F a s E 2 (x)\0), if this decay is known. The amplitude of this process is 
proportional to F, while it is cancelled in the ratio T(T(UU) -> hZ)/T{T{UU) -> 11). 
The conformal properties of QCD are especially important in transitions between lighter 
hadrons, e.g. (cc)- and (66)-bound states, where the multipole expansion is more sufficient 
[19]. In the case of (super)heavy quarkonium transitions, the Model I yields an enhancement 
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effect compared with the result coming from the Model II at ((a s /ir) G^„) ~ 0.012 GeV A . 
Such an enhancement of the decay width Y(T(UU) — > hZ) with guaranteed criterion c = 
(Ttot/^B) < 1 could have a quite significant effect on the shapes of bound-state resonances 
containing (super)heavy quarks in the search for new phenomena at the Tevatron and LHC. 
Unfortunately, the total decay width T tot is unknown. In Fig. 4, we show the dependence of 
the ratio of its total decay width T tot to the binding energy e#, i.e. c = T tot /e B on the real 
Yukawa couplings £ x q in the 2HDM at different values of x-boson mass. We observe that the 
binding critical ratio c is strongly sensitive to £ x q. Thus, we believe that in some kinematical 
regions there is a possibility for finding the fourth family quarks at the Tevatron's Run II (the 
LHC experiments will not miss it). If the scale of the Higgs-boson mass does not exceed 180 
GeV, the channels of Higgs decays to gg, 77 or 11 (p + p~ , r + r~) may give an enhancement 
effect which can be interpreted as the indication of heavy quark existence. For example, the 
decay of the CP-even lightest Higgs-boson h — > ##(77) with 114 GeV < rrih < 180 GeV 
could be enhanced (due to the 4th generation quark contribution) by a factor p ~ 8.9(4.5- 
5.5) for m 4 = 200 GeV and p ~ 8.5(4.0-5.0) for m 4 = 600 GeV. On the other hand, the decays 
of CP-even heavy Higgs H — > #(7(77) in the mass region 180 GeV < m H < 800 GeV could 
give 9 < p < 13(15 < p < 25) for m 4 =200 GeV and 9 < p < 27(15 < p < 57) for m 4 = 600 
GeV. This effect on the H Higgs-boson decay has a minor dependence of tan/5. We have 
shown that a reduction of a single [/-quark decay width leads to a significant enhancement 
of the signal of the T(UU) resonance. We investigated also the possible manifestation of the 
CP-even Higgs-boson states h and H in their rare decays h — > T(bb)^ and H — > T{bb)^ or 
even the production of more heavier bound states T(tt) in the rare decays of Higgs bosons 
H . We obtain that these decays can be detectable at the forthcoming experiments at the 
LHC. 

In the final state, the Higgs bosons h, the gauge bosons Z /W ± and even charged Higgs- 
bosons H ± are on mass-shell. Hence, the masses of heavy quarks can be reconstructed. As 
a trigger, one can choose a semileptonic decays U — > Dlui, U — > blvi, t — > blui (I = p,r) 
and the leptonic one like W — ■> Ivi with different lepton flavors to eliminate the backgrounds 
7*, Z, Z' — > //. On the other hand, we suppose that the efficiency for observing QQ events 
may be high enough because of the characteristic kinematics such as U — > DW, U — > bW 
and t — > bW. The comparison of the measured decay width mentioned above with theo- 
retically predicted ones can exclude or even confirm the fourth family fermions (quarks). 
Finally, related to our calculations here, some comments are order in the following: 

a. the (super)heavy quarks are considered to be non-relativistic; 

b. if the x Higgs-boson mass is m x ~ (9(100 GeV), then the decoupling limit (£ x q ~ C(l)) 
is appropriate only for uiq >2m t ; 

c. for the (tt)-bound state with £ x q ~ 0(1), the mass m x should be smaller than the lower 
bound given by the LEP 2 experiments [21]; 

d. if the heavy quarkonium mass is of the order 2mt, the x-boson contribution to the com- 
bined potential (4) becomes appreciable only for large values of t; xt > 6. 

We conclude that the (super)heavy quarkonia effects calculated here can be significant and 
should be considered seriously for searching for new physics beyond the SM. 
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